Creating Decoherence-Free Subspaces Using Strong and Fast Pulses 
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A decoherence-free subspace (DFS) isolates quantum information from deleterious environmental 
interactions. We give explícit sequences of strong and fast ("bang-bang", BB) pulses that create the 
conditions allowing for the existence of DFSs that support scalable, universal quantum computation. 
One such example is the creation of the conditions for collective decoherence, wherein all system 
partides are coupled in an identical manner to their environment. The BB pulses needed for this 
are generated using only the Heisenberg exchange interaction. In conjunction with previous results, 
this shows that Heisenberg exchange is all by itself an enabler of universal fault tolerant quantum 
computation on DFSs. 



Since the discovery of quantum error correcting codes 
(QECCs) jp, an arsenal of powerful methods has been de- 
veloped for overcoming the problem of decoherence that 
plagues quantum computers (QCs). A QECC is a closed- 
loop procedure, that involves freqüent error identifica- 
tion via non-destructive measurements, and concommi- 
tant recovery steps. Alternatively, decoherence-free sub- 
spaces (DFSs) and subsystems H, and dynamical 
decoupling, or "bang-bang" (BB) ]r|-pd|], are open-loop 
methods. A DFS is subspace of the system Hilbert space 
which is isolated, by virtue of a dynamical symmetry, 
from the system-bath interaction. The BB method is 
a close cousin of the spin-echo effect. All decoherence- 
reduction methods make assumptions about the system 
(S)-bath (B) coupling, embodied in a Hamiltonian of the 
general form H = Hs <S> Ib + Hsb + Is <& Hb- Here / is 
an identity operator and Hsb is the system-bath inter- 
action term, which can be expanded as a sum over linear, 
bilincar, and highcr order coupling terms: 
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Specializing to qubits, a typical assumption is p = 1, 

H i = 3 i -B i = °ï® B i, ( 2 ) 

ol— x,y,z 

where c?j = (of , of , of ) = (Xí,Y Í7 Zí) are the Pauli ma- 
trices acting on the i th qubit, and Bi = (Bf,Bf,B*) 
are arbitrary bath operators. This includes the indepen- 
dent errors model (all Bf different, a — x,y,z), and the 
collective decoherence model {Bf = B a Vi), important, 
respectively for the QECC and DFS methods. The bilin- 
ear terms are in general descibed by a second-rank tensor 
Gij , so that Hij = c?i ■ Gy • Sj . One of the main problems 
in applying the various decoherence-countering strategies 
is that, typically, the conditions under which they ap- 
ply are not wholly satisfied experimentally E ]. This 
problem is particularly severe for the DFS method, since 



it demands a high degree of symmetry in the system- 
bath interaction. Two main cases are known that admit 
scalable DFSs (i.e., subspaces that occupy a finite frac- 
tion of the system Hilbert space) : collective decoherence 
and the model of "múltiple qubit errors" (MQE) 
Collective decoherence, as defined above, assumes 
qubit-permutation-invariant system-bath coupling. This 
may be satisfied at ultralow temperatures in solid-state 
QC implementations, provided the dominant decoher- 
ence mechanism is due to coupling to a long-wavelength 
reservoir, e.g., phonons [G~4 15|. MQE assumes that the 
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system terms appearing in Hsb generate an Abelian 
group under multiplication (referred to below as the "er- 
ror group"). This is a somewhat artificial model that 
typically imposes severe restrictions on Bi and Gjj (ex- 
amples are given below). On the other hand, a two- 
dimensional DFS (encoding one logical qubit) can be con- 
structed using as few as 3 qubits under collective deco- 
herence conditions @, and unlike QECCs, requires no ac- 
tivo intcrvention other than the initial encoding and final 
decoding steps. One is thus faced with a rather frustrat- 
ing situation: the attractively simple DFS method im- 
poses symmetry demands that are likely to be perturbed 
in practice. Even though DFSs are robust with respect 
to such symmetry-breaking perturbations [pi|,|ïl| , and can 
be further stabilized by concatenation with QECCs |Ï7j , 
it is highly desirable to be able to artificially engineer 
conditions under which scalable DFSs may exist. Here 
we show how this can be accomplished for both the col- 
lective decoherence and MQE models, by combining the 
DFS encoding with the BB method. While such "envi- 
ronment engineering" methods have been proposed be- 
fore 0,^8), we show here how this can be accomplished 
for DFSs, assuming only physically reasonable resources. 
In particular, we show that by using decoupling pulses 
that are generated using only the isotropic Heisenberg 
exchange-interaction, one can transform the general lin- 
ear system-bath term Hi into a purely collective deco- 
herence term. Since the Heisenberg interaction is by itself 
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universal on the DFS-encoded qubits |[ïq , [L9| , this result 
has direct implications for the promising QC proposals 
that make use of the Heisenberg interaction to couple 
qubits 20 - 2|] , and in which the use of single-qubit oper- 
ations is preferably avoided ]2^-[25|]. 

Dynamical decoupling and DFS. — Let us start by 
briefLy reviewing the decoupling technique, as it pertains 
to our problem (for a thorough review see, e.g., p0[). 
Decoupling relies on the ability to apply strong and fast 
(BB) pulses ||, in a manner which effectively averages 
Hsb to zero. Since the pulses are strong one ignores 
the evolution under Hsb while the pulses are on, and 
since the pulses are fast one makes the short-time ap- 
proximation, i.e., exp[(A + B)t] « exp(Aí) exp(Bt) for 
[^4, £>] ^ 0. Systematic corrections are known M. The 
simplest example of eliminating an undesired unitary 
evolution U = cjíp(-üHsb) is the parity-kick sequence 
HI]] ■ Suppose we have at our disposal a fully controllable 
interaction generating a gate R such that "i? conjugates 
U" : RWR=W. Then the sequence UR)UR = I serves 
to eliminate U. 

Now, turning on the single-qubit Hamiltonian efXi 
for a time t = n/2ef generates the single-qubit gate 
Xi = i exp(— i^Xi). Each term in Hsb either commutes 
or anti-commutes with Xi since each term contains at 
least one factor of Xí,Yí,Zí. We call a term A "even" 
with respect to B if [A, B] = 0, "odd" if {A, B} = 0. If 
a term A in Hsb is odd with respect to (wrt) Xi then 
the evolution under it will be conjugated by the gate Xf. 
Xiexp(—iAr)Xi — exp(iAr). This allows for selectively 
removing this term using the parity-kick cycle, which we 
write as: [r, Xi, t, Xi\. Reading from right to left, this 
notation means: apply Xi pulse, free evolution for time 
r, repeat. Since every system factor in Hsb contains a 
single-qubit operator, it follows that we can selectively 
keep or remove each term in Hsb by using the parity- 
kick cycle. Note, however, that in general we have to 
use a short-time approximation since [U, R^UR] ^ @]. 
Further, without additional symmetry assumptions re- 
stricting p of Eq. this procedure, if used to eliminate 
all errors, requires a number of pulses that is exponential 
in N The reason is that without symmetry we will 
need at least two non-commuting single-qubit operators 
per qubit (e.g., X i} 

The DFS method uses a very different idea for over- 
coming decoherence Suppose that there exist 

states {IV'ï}} that are degenerate under the action of 
all system operators S a in the interaction Hamiltonian 
Hsb = J2 a s a ® B a , i.e.: S a \^i) = c a \tpi) Va, i, where 
c a are constants (more general conditions can be found 
||). Such a collection of states forms a subspace that 
aequires only an overall phase under the action of Hsb, 
and is therefore decoherence-free. The requisite degen- 
eracy arises from a symmetry in Hsb, such as collective 
decoherence. 

Symmetrization. — We now turn to showing how de- 



coupling may be used to create the conditions for DFSs. 
General group-theoretic arguments for using BB pulses 
for "symmetrizing" system-bath interactions, thus creat- 
ing DFS conditions, were given in pj. However, these 
proposals did not consider the MQE model and did not 
give an explicit Hamiltonian realization for the collec- 
tive decoherence model. Here we specialize to the MQE 
and collective decoherence models, and give explicit pulse 
sequences that respect the constraints imposed by phys- 
ically available resources. 

Generation of the MQE model. — The MQE model 
assumes that the system operators appearing in Hsb 
form an Abelian group Q under multiplication E3|. 

E.g., H S B = Eil^ z i z i+i ® B u or h sb = 
Et 7 ' E a=I , Sl! 0$í-iOm®B?. The dimension of the DFS 
supported by Q is 2 N /\Ç\, where \Q\ is the order of Q, 
which also counts the number of independent errors the 
DFS is immune to |13|. An Abelian group with M gen- 
erators has order 2- íí7 7Universal, fault-tolerant quantum 
computation can be performed on the MQE class of DFSs 
using the method of pq] , Briefly, this method involves 
a hybrid DFS-QECC approach, wherein lògic gates act- 
ing on DFS states are supplemented with fault-tolerant 
error detection and recovery. This active QECC inter- 
vention is needed since, unlike the collective decoherence 
case treated below, in the MQE case lògic gates take en- 
coded states on a trajectory that begins inside the DFS, 
leaves it, and then returns (as is also the case for compu- 
tation using QECCs). Let us now show how to generate 
the MQE conditions starting from a p — 2 system-bath 
Hamiltonian. 

We assume that single-qubit gates are available. In this 
case, it has been shown that the linear term, Hi, can 
be eliminated using 4 pulses, each acting simultaneously 
on all qubits |^|j9| . We reproduce this result and show fur- 
ther how single-qubit gates can efficiently decouple bilin- 
ear Hamiltonians Hij with nearest-neighbor interactions. 
Let H Dn = Y,í=i h í + and f7 nn = exp(-ÍTH nn ). 

Dcfine collective rotation operators 
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where N is even and R can be I, í or Z. First 
note that U nn = U nn (XU nn X) leaves only those lin- 
ear terms containing Xi, and all bilinear terms of the 
form afer™, YiZj. Let us apply Z to the outeome, i.e., 

U nn = U nn (ZU nn Z). This eliminates all linear terms 
in 4 pulses: [t,Y,t,X,t,Y,t,X] (where we have used 
Y = -íXZ). It also eliminates all Y%Zj <g> Bf-, leav- 
ing just ^i^i^j+l ® B f,j+i ( a = X ,V> Z )- This we 
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can rewrite as • f"cr" +1 

i.e., the even-numbered qubits act as baths for the odd- 
numbered qubits. Now let U nn = U nn (XoU nn Xo), which 
requires 8 pulses. At this point we are left only with er- 
rors of the form XjXj+\. These generate an Abelian 
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group denoted Q 2 x in |!6j. Since Q 2 x has N — 1 gen- 
erators its order is |Q2Jt| — '2 N ~ 1 , so that the DFS is 
2 N /2 N ~ 1 = 2 dimensional, i.e., supports a single en- 
coded qubit, which of course is not scalable. A larger 
number of encoded qubits can be supported by reduc- 
ing the dimension of the error group. This, in turn, 
rcquires a few more pulses. E.g., consider applying 
the BB pulse Z 3 Z^Z-jZ^ • • • to U nn . This climinatcs 
X 2 X 3 , XqXt, .... What is left, after 16 pulses, is 
the error group generated by {^j-i^il^i! denoted 
Qx hi [JL3|. It has order \Qx\ — 2 N ^ 2 , thus supporting 
a 2 N /2^ = 2 N ' 2 dimensional DFS. This DFS encodes 
N/2 qubits, so it is scalable. The methods of p(| now 
apply for the purpose of fault-tolerant universal quantum 
computation using the hybrid DFS-QECC method. 

Note that we can also go further and eliminate all 
second order coupling terms: U nn (ZoU nn Zo) — I and 
also uses a total of 16 collective pulses. If there 
is a next-nearest-neighbor system-bath interaction, it 
can similarly be removed using the collective operator 
Roo — R1R5R9R13 ■ ■ ', etc. for longer-range interac- 
tions. These manipulations will leave higher order MQE 
models; which option to choose will depend on which 
pulse sequences are most easily implementable. It should 
be clear that this method of generating MQE models is 
quite general: given a system-bath Hamiltonian, one can 
design a set of BB pulses that will transform this Hamil- 
tonian into a desired Abelian error group. The number 
of pulses will scale with the system-bath coupling order 
p [Eq. (0)] and the interaction range r (r = 1 for nearest- 
neighbors, etc), and we have given explicit sequences for 
the case p, r < 2. We note that this combination of de- 
coupling with a hybrid DFS-QECC strategy is, as far as 
we know, the first time that all three methods for com- 
batting decoherence have been used together. 

As a final comment on generating MQE models, we 
note that the analysis above applies also to the case 
where it may be preferable to control two-qubit "prod- 
uct" Hamiltonians of the form XiXj and Y{Yj. Simi- 
lar to the case of controllable single-qubit gates, we now 
have the gates X^Xj = iexp(—ÍTrXiXj/2) and YiYj — 
i exp(—ÍTrYiYj/2). Such gates could be implemented nat- 
urally, e.g., in certain superconducting QC implementa- 
tions j27|. It is simple to check that the product Hamil- 
tonians can be used to decouple any linear system-bath 
Hamiltonian, and any bilinear term other than af oj*. In 
fact, we can construct the R gates [Eq. (||)] by simply 
turning on all nearest neighbor gates cr"cx" +1 , and the 
Ro gates by simultaneously turning on all next-nearest 
neighbor gates o-"ct" +2 . Therefore the methods used 
above for the case of single-qubit gates apply directly. 

Generating collective decoherence. — We now turn to 
collective decoherence. To do so, we consider the impor- 
tant case of a controllable Heisenberg exchange Hamilto- 
nian JijSi ■ (Tj, crucial for the operation of, e.g., quantum 



dot QCs [£0| an d donor atom nuclear (2ÏJ or electron [£2| 
spin QCs. We will show that using a few collective pulses 
generated by the Heisenberg interaction alone, we can 
symmetrize any linear system-bath interaction $i ■ Bi, 
such that only a block- collective component remains. This 
collective decoherence can then be avoided by encoding 
into a 4-qubit DFS or a 3-qubit DFS/noiseless sub- 
system (NS) f§. Let 

Oij = exp(— ínài ■ cr.,/4) 

A simple calculation shows that Oij is a SWAP gate 
for the Pauli matrices: O^SiOij ~ tjj. From this follows 

0\j (<Ji ± (Tj) Oij = ± (Si ± <jj), i.e., all differences Oi — dj 
are odd wrt Oij, and hence can be eliminated. Let us 

then rewrite °i ' = J2p=± Z)j=i (^2j + P<?2j-l) ' 



B$j, where B± 



(B 2j ± B 2 j-i)/2. To elimi- 
nate the nearest-neighbor differences (c?2j — <?2j-i) we 
can use the collective BB pulse O — ^2j-i,2j- 
This leaves only the (block-)collective decoherence term 
Y]j—i (<?2j + (?2j-i) • Bj~, which in turn we can rewrite 

aS E/3=± E^íl 2 í^2j+2 + <?2j + l + (<?2j + 0*23-1)) ■ Btfi 



where B. 



+,± 
23 



(B~2j + 2 i B£j)/2. We can now eliminate 



the next-nearest neighbor differences (í?2j+2 — &2j) and 
(0*23+1 — oVz-i) using a second collective pulse Oo — 
®j=i _1 023-1,23+1 023,2i+2· At this point we are left 
just with the collective decoherence terms on blocks of 
4 qubits, and the encoding into the 4-qubit DFS be- 
comes relevant. This scheme uses a total of 6 collective 
BB pulses: [0 ,0,t,0^ ^,0^,0, t^O^]. The apparent 
drawback of needing next-nearest neighbor interactions 
(in the Oo pulse) can be avoided by swapping local gates, 
e.g.: Vi+2 = Ol +líi+2 O iíi+í O i+li i +2 , at the expense of 
more pulses. Note that in a 2D hexagonal arrangement 
i, i + 1 and i, i + 2 can all be nearest neighbors. 

We can also symmetrize into blocks of 3, which 
can be used for the 3-qubit DFS/NS |]. Let us 
rewrite H 3 = J^L·i <?i ' B i = (o'l + 0*2 + o* 3 ) ■ À + + 
(5*2 — 0*1) • À~ + 5*3 • C, where À ± = (B 2 ± Bi)/2 
and C = B3 — A + . We can eliminate 0*2 — 0*1 us- 
ing Oi 2 : Ui(t) = oíp(-íH 3 t)0\ 2 exp(-iH 3 T)Oi 2 
1 x [ ) — 2ít(((t 1 +0*2 + 0*3) • A + + 0*3 • C . Next consider 

U 2 (r) = U 1 (T/2)Ol 3 U 1 (r)0 23 

— e -iT[3{3 1 +3 2 +3 3 )-À + + (23 2 +8 3 )-C] 

Finally, 

U 2 (t)0\ 2 U 2 (t)0 12 = e -^( í i+^+^)-(Bi+s 2 +s 3 ) i 

leaving only the collective component. This scheme uses 
a total of 14 collective BB pulses. 

To quantum compute universally on DFSs it is nec- 
essary to couple blocks of DFS qubits, in order to 
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implement a controlled- lògic gate |ÏJ,|Ï9|]. An ex- 
tra symmetrization step is thus required, creating col- 
lective decoherence conditions over blocks of 6 or 
8 qubits. This is a simple extension of the pro- 
cedure above. E.g., for two coupled 4-qubit DFS 
blocks, we need collective pulses of the form Ooo = 

Ç&f=l 3 02j-l,2j+302j,2j+i02j + l,2j+502j+2,2j+6- By 

swapping local gates we can again avoid direct control 
over long-range interactions. The corresponding increase 
in the number of gates may well be a worthwhile trade- 
off. Similar pulse sequences can be found for creating 
block-collcctivc decoherence conditions over 6 qubits, for 
computation using the 3-qubit NS. 

Discussion and conclusions. — The prospect of 
decoherence-free quantum computation is very attrac- 
tive, but so far ideas for obtaining the conditions en- 
abling the existence of scalable decoherence-free sub- 
spaces (DFSs) focused mostly on lowering the tempera- 
ture and neglecting other sources of decoherence [14 15 



The exception are previous existential results showing 
how decoupling methods can be used for symmetrization 
of system-bath interactions JÏ0| . In this work we showed 
explicitly how conditions for the two most important ex- 
amples of DFSs (the models of collective decoherence 
and múltiple qubit errors ]Ï3| ], MQE) can be actively 
generated using symmetrizing cycles of fast and strong 
decoupling ( "bang-bang" , BB) pulses. In the MQE case 
a cycle of 16 pulses suffices to symmetrize a system-bath 
Hamiltonian with arbitrary linear and nearest-neighbor 
bilinear couplings, such that conditions enabling the ex- 
istence of a scalable DFS are established. This result is 
applicable for quantum computer (QC) proposals where 
single-qubit gates are easily tunable. In this case a hybrid 
DFS-active quantum error correction scheme developed 
in p6f , can be used for universal, fault-tolerant quantum 
computation. 

In the case of collective decoherence a very attrac- 
tive general picture is emerging, from a combination of 
previous studies and this work. The collective deco- 
herence model was first proposed as an example allow- 
ing the existence of DFSs with the property of a scal- 
able encoding It was later realized that universal 
quantum computation is possible on these DFSs using 
the Heisenberg exchange interaction alone p5| , ^9[ . Our 
present result shows how collective decoherence condi- 
tions can be actively created using a few pulses generated 
using only Heisenberg exchange. Since our method relies 
on BB pulses some degree of leakage out of the DFS 
(due to imperfect symmetrization) is inevitable. For- 
tunately, such leakage errors can also be reduced using 
Heisenberg-generated BB pulses PH , or detected with a 
circuit that utilizes, again, only Heisenberg exchange M. 



The combination of all these results shows that Heisen- 
berg exchange is all by itself an enabler of universal fault 
tolerant quantum computation on decoherence-free sub- 
spaces. This has potentially important applications for 
those solid-state proposals of QCs where Heisenberg ex- 
change is the natural qubit-qubit coupling mechanism 



p] 22 1. The combination of the decoupling method with 
encoding methods developed in the quest to protect frag- 
ile quantum information thus seems to be a promising 
route towards robust implementations of QCs. 
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